A theory of amorphous packings of binary mixtures of hard spheres 
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We extend our theory of amorphous packings of hard spheres to binary mixtures and more gener- 
ally to multicomponent systems. The theory is based on the assumption that amorphous packings 
produced by typical experimental or numerical protocols can be identified with the infinite pressure 
limit of long lived metastable glassy states. We test this assumption against numerical and experi- 
mental data and show that the theory correctly reproduces the variation with mixture composition 
of structural observables, such as the total packing fraction and the partial coordination numbers. 



Amorphous packings of hard spheres are ubiquitous 
in physics: they have been used as models for liquids, 
glasses, colloidal systems, granular systems, and pow- 
ders. They are also related to important problems in 
mathematics and information theory, such as digitaliza- 
tion of signals, error correcting codes, and optimization 
problems. Moreover, the structure and density (or poros- 
ity) of amorphous multicomponent packings is important 
in many branches of science and technology, ranging from 
oil extraction to storage of grains in silos. 

Despite being empirically studied since at least sixty 
years, amorphous packings still lack a precise mathemati- 
cal definition, due to the intrinsic difficulty of quantifying 
"randomness" [1]. Indeed, even if a sphere packing is a 
purely geometrical object, in practice dense amorphous 
packings always result from rather complicated dynam- 
ical protocols: for instance, spheres can be thrown at 
random in a box that is subsequently shaken to achieve 
compactification [2] , or they can be deposited onto a ran- 
dom seed cluster [3]. In numerical simulations, one starts 
from a random distribution of small spheres and inflates 
them until a jammed state is reached [4, 5]; alternatively, 
one starts from large overlapping spheres and reduces 
the diameter in order to eliminate the overlaps [6-9]. In 
principle, each of these dynamical prescriptions produces 
an ensemble of final packings that depends on the de- 
tails of the procedure used. Still, very remarkably, if the 
presence of crystalline regions is avoided, the structural 
properties of amorphous packings turn out to be very 
similar. This observation led to the proposal that "typi- 
cal" amorphous packings should have common structure 
and density; the latter has been denoted Random Close 
Packing (RCP) density. The definition of RCP has been 
intensively debated in the last few years, in connection 
with the progresses of numerical simulations [1, 10]. 

Nevertheless, the empirical evidence, that amorphous 
packings produced according to very different protocols 
have common structural properties, is striking and call 
for an explanation. This is all the more true for binary or 
multicomponent mixtures, where in addition to the usual 
structural observables, such as the structure factor, one 
can investigate other quantities such as the coordination 
between spheres of different type, and study their varia- 



tion with the composition of the mixture. 

In earlier attempts to build statistical models of pack- 
ings, only the main geometrical factors, such as the rel- 
ative size and abundance of the different components, 
were taken into account [11-13]. More precisely, these 
models focus on a random sphere in the packing and 
its first neighbors, completely neglecting spatial correla- 
tions beside the first shell and all the global geometric 
constraints. This already accounts for the main qualita- 
tive structural properties of random packings. However, 
in order to obtain a quantitative description, some free 
parameters have to be introduced and adjusted to match 
with experimental data. 

To go beyond these simple models, many authors pro- 
posed that random packings of hard spheres can be 
thought as the infinite pressure limit of hard sphere 
glasses [14-20]. This is very intuitive since a glass is 
a solid state in which particles vibrate around amor- 
phous reference positions, and vibrations are reduced on 
increasing pressure. A typical algorithm attempting to 
create a random packing starts at low density and com- 
presses the system at a given rate. During this evolu- 
tion, when the density is high enough, relaxation becomes 
more and more difficult until at some point the system 
is stuck into a glass state [19, 20]; at this point further 
compression will only reduce the amplitude of the vibra- 
tions. In a nutshell, this is why amorphous packings can 
be identified with glasses at infinite pressure. 

The main advantage of this identification, if it holds, is 
that a glass is a metastable state that has a very long life- 
time; therefore, its properties can be studied using con- 
cepts of equilibrium statistical mechanics. In this way a 
complicated dynamical problem (solving the equations of 
motion for a given protocol) is reduced to a much more 
simple equilibrium problem. In [18, 20] it was shown, in 
the case of monodisperse packings, that this strategy is 
very effective since it allows to compute structural prop- 
erties of random packings directly from the Hamiltonian 
of the system, without free parameters and in a controlled 
statistical mechanics framework. Note that the existence 
of an equilibrium glass transition in hard sphere systems 
(or in other words the existence of glasses with infinite life 
time) has been questioned [5, 21]. Although very inter- 
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esting, this problem is not relevant for the present discus- 
sion since we are only interested in long-lived metastable 
glasses that trap dynamical algorithms. At present it 
is very well established by numerical simulations [22, 23] 
that for system sizes of N < 10 4 particles and on the time 
scales of typical algorithms, metastable glassy states ex- 
ist, at least in d > 3. This is enough to compare with 
most of the currently available numerical and experimen- 
tal data. Finally, the relation of this approach to special 
packings such as the MRJ state [1] and the J-point [8] 
has been discussed in detail in [20]. 

The aim of this paper is to extend the theory of [18, 20] 
to binary mixtures. This allows to compare quantita- 
tively the predictions of the theory and the results of 
numerical simulations. We will focus in particular on the 
variation of density and local connectivity as a function 
of mixture composition. These results constitute, in our 
opinion, a stringent test of the assumption that random 
packings reached by standard algorithms can be identi- 
fied with infinite pressure metastable glasses. 

Theory - The equilibrium statistical mechanics com- 
putation of the properties of the glass is based on stan- 
dard liquid theory [24] and on the replica method [17, 25] 
that has been developed in the context of spin glass the- 
ory [26, 27]. For monodisperse hard spheres, it has been 
described in great detail in [20] . The extension to multi- 
component systems is straightforward following [28]; de- 
tails are given in Appendix. 

Here we just recall some features of this approach, 
based on simple physical considerations. The basic as- 
sumption of the method are that i) crystallization and 
phase separation are strongly suppressed by kinetic ef- 
fects so that the liquid can be safely followed at high den- 
sity, and it) that at sufficiently high density, the liquid is 
a superposition of a collection of amorphous metastable 
states. Namely, in the liquid, the system spends some 
time inside one of these states, and sometimes under- 
goes a rearrangement that leads to a different state [29]. 
Each state is characterized by its vibrational entropy per 
particle, denoted by s, and the number of such states is 
assumed to be exponential in N , so defining a conjigura- 
tional entropy E(y>, s) — N~ 1 logJ\f(ip, s), being N(<p, s) 
the number of states having entropy s at density p. On 
increasing the density, the liquid is trapped for longer and 
longer times into a metastable state, until at some point 
the transition time becomes so long that for all practical 
purposes the system is stuck into one state: it then be- 
comes a glass. To compute the properties of the glassy 
states, the central problem is to compute the function 
E(y>, s). This can be done by means of a simple replica 
method introduced by Monasson [27] . One introduces m 
copies of each particle, constrained to be close enough, 
in such a way that they must be in the same metastable 
state. Then, the total entropy of the system of m copies 
is given by <S(m, ip\s) — ms + Yi(ip, s); the first term gives 
the entropy of the m copies in a state of entropy s, while 
the second term is due to the multiplicity of possible 
states. The total entropy of the system at fixed density 



0.75 




0.65 



0.2 0.4 0.6 0.8 

Volume fraction of small component 



FIG. 1: (Color online) Packing fraction <pj as a function of 
77 = 1/(1 + xr 3 ) at fixed r. Full symbols are numerical data 
from this work. Open symbols are experimental results from 
Ref. [30]. Lines are predictions from theory, obtained fixing 
Ej = 1.7. Note that the large r-small 77 region could not be 
explored, since for such very asymmetric mixtures the large 
spheres form a rigid structure while small spheres are able to 
move through the pores and are not jammed [12, 13]. 



tp is obtained by maximizing over s, i.e. 

S(m, ip) — max[E(</9, s) + ms] = £(y>, s*) + ms* , (1) 

s 

where s*(m,p) is determined by the condition 
S S(<^, s) = m. Then it immediate to show that 



s*(m, tp) = d m S(m, tp) , 

E(y>, s*(m, tp)) — ~m 2 d m [S(m, tp)/m] 



(2) 



The knowledge of S(m, ip) allows to reconstruct the curve 
E(y?, s) for a given density by a parametric plot of Eqs. (2) 
by varying m. The function E(y>, s) gives access to the 
internal entropy and the number of metastable glassy 
states; from this one can compute their equation of state, 
i.e. the pressure as a function of the density; in particu- 
lar, for each set of glassy states of given configurational 
entropy £j, one can compute the density ipj (jamming 
density) at which their pressure diverges. Since tpj turns 
out to depend (slightly) on Ej, a prediction of the the- 
ory is that different glasses will jam at different density: 
amorphous packings can be found in a finite (but small) 
interval of density [19, 20]. 

Results for binary mixtures - The details of the compu- 
tation of the function S(m, ip) for a general multicompo- 
nent mixture, based on [20], can be found in Appendix. 
Here we consider a binary mixture of two types of spheres 
fi = A, B in a volume V, with different diameter and 
density p M = N^/V. We define r = D A /D B > 1 the 
diameter ratio and x = Na/Nb the concentration ratio; 
Vs(D) = ttD 3 /6 the volume of a three dimensional sphere 
of diameter D; ip = paV^(Da) + PbV?,(Db) the packing 
fraction; r\ = pbVj,(Db)/p — 1/(1 + xr 3 ) the volume 
fraction of the small (B) component. 
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FIG. 2: (Color online) Partial average coordination numbers (small-small, small-large, large-small, large-large) as a function of 
volume fraction of the small particles n — 1/(1 + xr 3 ) for different values of r. Full symbols are numerical data from this work. 
Open symbols are experimental data from Ref. [31]. Note that in the lower right panel a different scale is used for zi 3 . 



Once an equation of state for the liquid has been cho- 
sen, the jamming packing fraction ipj is given in term 
of Ej by the solution of Ej(<^) = Ej. The average co- 
ordination numbers at tpj are denoted z^fyj), but we 
checked that their variations with ipj are negligible. We 
used in d = 3 the equation of state proposed in [32], us- 
ing the Carnahan-Starling equation for the monodisperse 
system [24]. The latter, as well as £j(<p) and Zn V (<pj), are 
given respectively in Eqs. (20), (16), (19), in Appendix. 

Numerical simulations - We produced jammed pack- 
ings of binary mixtures of N = 1000 hard spheres using 
the code developed by Donev et al. [33, 34]. In this al- 
gorithm spheres are compressed uniformly by increasing 
their diameter at a rate dD/dt = 2j while event-driven 
molecular dynamics is performed at the same time. In or- 
der to obtain a perfectly jammed final packing, the later 
stages of compression must be performed very slowly. 
On the other hand, at low density slow compression is 
a waste of time, since the dynamics of the system is very 
fast. Following [35], we find a good compromise by per- 
forming a four stages compression: starting from random 
configurations at tp = 0.1, i) the first stage is a relatively 
fast compression (7 = 10 -2 ) up to a reduced pressure 
p = (3P/p = 10 2 ; then we compress at ii) 7 = 10~ 3 up 
to p = 10 3 ; in) 7 = 10~ 4 up to p = 10 9 ; w) 7 = 10~ 5 
up to p = 10 12 . The first stage terminates at a density 
(p ~ 0.6, and is fast enough to avoid crystallization and 
phase separation. During the following stages the system 
is already dense enough to stay close to the amorphous 
structure reached during the first stage. Still, little rear- 



rangements (involving many particles) are possible and 
allow to reach a collectively jammed final state [35]. In 
the final configurations, we observe a huge gap between 
contacting (typical gap ~ 10~ n Z?) and non-contacting 
(typical gap > 10~ 6 £>) particles. We then say that two 
particles are in contact whenever the gap is smaller than 
10~ 8 £>. Typically, a small fraction (< 5% of the total) 
of rattlers, i.e. particles having less then 4 contacts, is 
present. Once these are removed, the configuration is 
isostatic (the total number of contacts is 6N) within 1% 
accuracy. 

Comparison of theory and numerical/ experimental 
data - During the four stages of compression, the pressure 
initially follows the liquid equation of state up to some 
density close to the glass transition density [20, 36, 37]. 
Above this density, pressure increases faster and diverges 
on approaching jamming at ifj. The exact point where 
this happens depends on compression rate. This is a nice 
confirmation of a prediction of the theory, that differ- 
ent glassy states jam at different density; it was already 
observed in [5, 35] and recently discussed in great de- 
tail in [23, 38]. Within the theory ifj is related to T,j, 
the value of configurational entropy at which the system 
falls out of equilibrium; hence there is one free parame- 
ter, Ej, that depends on the compression protocol. The 
equation of state of the glass obtained numerically with 
our protocol corresponds within our theory to Ej ~ 1.5. 
We decided to use the value Ej = 1.7 that gives the best 
fit to the numerical data. This is consistent with previ- 
ous observations, that the configurational entropy is close 
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to 1 when the system falls out of equilibrium [39] . A de- 
tailed discussion of the behavior of pressure can be found 
in Appendix. 

In figure 1, we report the jamming density for differ- 
ent mixtures, putting together our numerical results and 
experimental data from Rcf . [30] , and the theoretical re- 
sults. Note that a single "fitting" parameter Ej, that is 
strongly constrained, allows to describe different sets of 
independent numerical and experimental data. The pre- 
diction of our theory are qualitatively similar to previous 
ones [12, 13], but the quantitative agreement is much bet- 
ter. Interestingly, a similar qualitative behavior for the 
glass transition density has been predicted in [36, 37]; 
although there is no a priori reason why the jamming 
and glass transition density should be related [23], it is 
reasonable to expect that they show similar trends [37] . 

Finally, in figure 2 we report the average partial con- 
tact numbers for different mixtures. These values have 
been obtained by removing the rattlers from the pack- 
ing. As discussed above, the total coordination is close 
to the isostatic value z = 6. This is also a non trivial 
prediction of the theory, see Appendix. As it can be seen 
from figure 2, the computed values agree very well with 
the outcome of the numerical simulation, at least for r 
not too large. Some discrepancies are observed in the 
contacts of the large particles for large r. We produced 
packings of N = 10 4 and checked that these are not finite 
size effects. Also, inspection of the configurations seem 
to exclude the presence of phase separation. However, 
for these values of r and x, a large fraction of rattlers 
(~ 10%) is present within the small particles. This might 
affect the determination of the partial contacts. It would 
be interesting to check if better results are obtained us- 



ing different algorithms. Experimental data from [31] are 
also reported in the right panel of figure 2. 

Conclusions - In this paper we have extended our the- 
ory of amorphous packings to binary mixtures and have 
tested it against numerical and experimental data. In 
particular we have shown that the theory correctly pre- 
dicts the variation of total density (or porosity) and lo- 
cal coordination with mixture composition. We have also 
shown that the behavior of pressure during compression 
follows the predictions of the theory. A striking predic- 
tion of the theory is that different compression proce- 
dures lead to different final densities, which seems to be 
confirmed by numerical data, see also [23, 38]. Note that 
the only free parameter is the value of Ej, that is still 
strongly constrained (it must be close to 1). It affects 
slightly the values of density (by varying Ej in the rea- 
sonable range one can change ipj of ~ 10%, see figure 3 in 
Appendix) and does not affect at all the curves in figure 2 
for the local coordinations. As stated in the introduction, 
we believe that these results constitute a stringent test 
of the idea that amorphous packings can be considered 
as the infinite pressure limit of metastable glassy states. 
Note that our results have no implications on the exis- 
tence of an ideal glass transition [5, 21]. Indeed, we only 
used numerical data obtained using fast compressions. It 
is possible that using much slower compressions physics 
changes dramatically and the glass transition is avoided. 
Still, these time scales are out of reach of current algo- 
rithms and experimental protocols. 
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Appendix 



I. BEHAVIOR OF THE PRESSURE DURING COMPRESSION 
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FIG. 3: (Color online) Inverse reduced pressure, p/(J3P) as a function of the packing fraction ip for a mixture with r = 1.4 and 
x — 1. Numerical data are obtained using two different protocols. In the first, compression is started at low density. In the 
second, compression is started from an equilibrated configuration at ip = 0.58. The equation of state of different metastable 
glasses, corresponding to different Ej, are reported as dashed lines. The dot-dashed line is the pressure of the ideal glass, 
corresponding to = 0. A numerical estimate [23] of the Kauzmann pressure, pk = 34.4, is reported as a dotted horizontal 
line. 



Here we discuss in detail the behavior of pressure during compression in our numerical simulations. In figure 3 we 
report the evolution of the inverse reduced pressure during the first and second stages of compression {Numerical 1) 
for a mixture with x = 1 and r — 1.4. We observe, as in similar studies [5, 23, 35, 38], that the pressure follows the 
liquid equation of state up to some density, that for this system is around ip ~ 0.56. At this value of density, the 
relaxation time of the liquid becomes long enough that liquid relaxation is effectively frozen on the compression time 
scale used in this work and the system falls out of equilibrium. Although we did not measure relaxation time directly, 
this has been done in related works that confirmed the correctness of this statement [23, 37, 38]. Above this density, 
pressure increases faster and diverges on approaching jamming around <pj ~ 0.65. The numerical equation of state is 
compared with that of a glass state corresponding to Sj = 1.5. This is consistent with previous observations, that 
the configurational entropy is close to 1 when the system falls out of equilibrium [39]. In the same plot, we report 
the curve (Numerical 2) obtained starting the first stage from a carefully equilibrated liquid configuration of the same 
mixture at <p = 0.58, kindly provided by L. Berthier (see [23] for details on how this configuration was produced and 
equilibration was checked). In this case, since the relaxation time of the liquid at that density is already very long 
compared to our compression rate, the system falls immediately out of equilibrium and the pressure increases fast until 
jamming occurs at a higher density compared to the previous case. In our interpretation, this corresponds to a glassy 
state with lower Ej (compare with the theoretical curve for Ej ~ 1.2). This is a nice confirmation of a prediction of the 
theory, that different glassy states jam at different density. Finally, we report, for the same system, the numerically 
extrapolated value of the ideal glass transition pressure, px — 34.4, see Ref. [23] for details of the careful extrapolation 
procedure. Again, this corresponds well (within 10%) to the computed value px = 31.8 from the theory, see figure 3. 
Note that this coincidence does not prove the existence of the Kauzmann transition, since on much larger time and 
length scales than the ones explored in these simulations, a crossover to a non mean- field behavior might happen. 
Still, the coincidence shows that, on the length and time scales explored by current numerical simulations, the full 
mean field phenomenology is observed, including the apparent extrapolation of Kauzmann transition. Whether this 
transition really exists remains a major open point of the field, that however does not affect the results presented 
here. 
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II. REPLICA THEORY FOR MULTICOMPONENT MIXTURES 



Here we show in detail how to generalize the computation of [20] to multicomponent mixtures. We will not repeat 
the discussion of [20] but we only explain how to modify it for mixtures. Reading Appendix B and C and section VII 
of [20] is necessary to follow the discussion. 

We consider a multicomponent system in a (large) volume V with partial densities p M = N^/N. The total density 
is p = an d we define x^ = Py.jp. Spheres of type p have diameter Z? M and we define = (D^ + D v )/2. 

We denote by Ct d = 2ir d / 2 /T(d/2) the d-dimensional solid angle, and by V d (D) = ir d / 2 D d /T(l + d/2) the volume of a 
d-dimensional sphere of diameter D. The hard sphere potential ^^(r) is infinite for r < and zero otherwise: we 
denote Xnv{ r ) = ex P[~ ( ?W( r )]- 5/^( r )j as usual, is the p-v pair correlation function [24]. We will use the shorthand 
notations g^ v = g^D^) for the contact values of <7 M1 ,(r), and V^ v = Vd{D lxv ) for the volume of a sphere of diameter 
/>,,,- 

As discussed in [28], we assume that in the replicated liquid molecules are built of particles of the same type. This 
amounts to assume that in a glassy state particles of different type cannot easily exchange. However, as the diffusion 
constant is always finite in a glass, this assumption is wrong, since particles can exchange also within a state. This 
is all the more true in the liquid phase before the glass transition. Therefore, in order to correct this error, we will 
subtract the mixing entropy from the configurational entropy we will compute in the following, in order to obtain the 
correct physical result. 

According to [28], to each molecule of the replicated liquid we can attach a label p = 1, ■ ■ ■ , n according to the 
type of particles that build the molecule. We denote coordinates in a molecule by x = (x\, ■ ■ ■ , x m ). We have then to 
compute the free energy functional of a liquid made of n species of molecules, with a single molecule density PlJ ,(x) 
and an interaction x^(x, y) = X^( x a - Da)- 



Effective liquid 



We assume that the vibrations of the m copies are described by a Gaussian distribution, corresponding to harmonic 
vibrations: 

(x n -X) 2 

/m 2A 

The widths are variational parameters and we will maximize the entropy with respect to them at the end. We 
choose a replica (say replica 1) as a reference and consider the vibrations of the other m — 1 particles around the 
reference one. We want to construct an expansion assuming that is small. For A^ = 0, the m — 1 copies coincide 
with the reference one, and S(m, ip) is given by the entropy S(tp) of the non-replicated liquid plus the entropy of m — 1 
harmonic oscillators of spring constant A^ (averaged over the concentration of the different species) : 



5 (0) (m, <p) = S(ip) + x»S harm {m, A„) . (4) 

The harmonic part of the entropy is computed straightforwardly from the term X^u / dxpn(x)[l — logp /i (a;)] in the 
free energy functional, see section V of [20], where the function 5/j arm (m, A) is also defined. 

A first order approximation is obtained by considering the effective two-body interaction induced on the particles of 
replica 1 by the coupling to the m—l copies. This can be justified on the basis of a diagrammatic expansion following 
the derivation in Appendix B of [20], see also [15]. It is possible to show [24] that the diagrammatic expansion of 
the free energy functional for a multicomponent liquid is the same as the one of a simple liquid, provided a label pi 
is attached to each vertex i, and in such a labeled diagram p^Xi) is placed on each vertex, x^y. ( x i ~ x j) on each 
link, and a sum over pi is performed in addition to the integration over X{. The same is true for the molecular liquid. 
Hence, the treatment of diagrams performed in Appendix B of [20] can be repeated exactly for a mixture, taking 
into account the presence of the additional labels. This is straightforward and leads to the introduction of effective 
potentials that now depend on the type of particles involved. In particular, one obtains: 

/m 
dx 2 ■ ■ ■ dx m dy 2 ■ ■ ■ dy m p^x) Pv {y) J[ e -*M*(*.-¥.) 
o=i (5) 

= e -4>M*i-vi)[ 1 + Q^ Xl - yi )} . 
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A first order approximation to S(m, tp) is then obtained by substituting to the entropy S(ip) of the simple hard 
sphere mixture the free energy of a liquid of particles interacting via the potential 4> e JJ (r) : 

S^(rn, i p) = -F[ V ,cl> e /J(r)}+J2^Sharrn(rn,A fl ) . (6) 

fj- 

It is now evident that we can obtain better approximations of the true function S(m, tp) by considering also the three 
body interactions induced on particles of the replica 1, and so on. Here, for simplicity, we will limit ourselves to 
consider the two-body interactions. To the same approximation, the pair correlation function g^ v (r) of the system is 
given by the correlation of the effective liquid of particles interacting via </> M „(r). For small A^, (j) e JJ is close to M „ 
and the free energy in Eq. (6) can be computed in perturbation theory around the normal liquid [24], which is taken 
as an input for the calculation of the properties of the glass. 



Small cage expansion 



Now we follow the discussion in section VII of [20] to compute the entropy of the effective liquid. We assume that 
the cage radius is small, therefore can be considered as a perturbation, and one has (jfJJ \r) <~ (f>^ u {r) — Q^{r). 
Using the general relation [24] 

SFjtp^^jr)} _ PllPv 

we obtain 



-F[<p, 4> e //(r)} ~ S(<p) + J2 ^ J dr 9^ 



(r)Q^(r) , (8) 



where S(<p) = —F{tp, ^^(r)] is the entropy of the liquid. The function Q^ u {r) is different from zero only for r — D^ u <~ 
0(y/A^ + A v ). If yJA^ + A v <§; D^ v we can assume that the function g^v^f) is essentially constant on the scale 
yJA^ + A v , then g^{r) ~ (r) for r- ~ 0( y /A, 1 + A v ), and 

J drg^^Q^r) ~ g^ j drx^{r)Q^{r) . (9) 

To compute the previous expression, all the considerations of Appendix C of [20] can be repeated. Instead of a 
single function Q(r), one has Q^ v {r), and in the two Gaussians one has A^ and A v instead of A. At the very beginning 
of the discussion, then, one finds that it is enough to replace 2A — > A^ + A v and follow the same steps. We obtain 

J drx„„(r)Q„„(r) = _i_^ 2 (^ + A v ) Q (m) . (10) 

where Qo{m) is defined in Appendix C of [20]. 
Putting all together, we find the final result 

S(m, tp- {Afj}) = x^Sharmim, + S(ip) + ^ ^g^-^J2(Ay + A v ) Q (m) (n) 



This expression generalizes the result of section VII of [20] to multicomponent mixtures; it has to be optimized over 
all the A^ to get the replicated entropy S(m, tp). 



Replicated free energy for binary mixtures 



The optimization with respect to A^ leads to n coupled equations that are not easy to solve in general. Therefore 
in the following we focus on the case of binary mixtures that is of interest here. 
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For convenience we denote A a = A and As = B the cage radii of the two species. We then obtain 

dS{m,(p;A,B) _ p A g AA V i f A dQo(m) PApBgABV^d Q (m) d(l - m) _ 

OA ~ 2 P D AA s/A pD AB ^/2(A + B) ^ 2A 

dS{m,y-A, B) = p 2 B 9BBV d BB dQ {m) p A p B 9ABV d AB d Q {m) _ ^ d(l — m) _ 

8B 2pD BB ^B pD AB ^2(A + B) XB 2B 



(12) 



xr^ 9 aa - V5g BB + V2 1 ~ g f i±I ) 5AB =0, (13) 



Defining 5 = B/A, x = pa/pb and r = Da/D b , we obtain the following equation for 5: 

) fi+ry- 1 

and from its solution we obtain the optimal values of A and B = 5 A: 

1-m (£ d )(l + .x) 



2Q (m) 2>T 



i ^7/ > ( • 9AB - i9 - Ui 



(14) 



' ^^2(1 + 5) 



where the packing fraction is <p = paV<i{Da) + pBVd{D B ) and (D p ) = x^D?. Substituting these in S(m, tp; A, B), 
we finally get 

a \ a i \ d. d. d 1 ... 2 d d<p Qo(m)\/A*A 

S(m, = Sfo>) - -(1 - m - logm) - -(1 - m) log(27rA*) - ^ jj^ 1 ~ m ) lo sW + 7^y (1 + ^2 (15) 

A = x 2 g AA D d AA 1 + xg AB D d A! 1 ^2{l + 5) + g BB D d B ^5 



From this expression one can compute the complexity using Eq. (2) [20]; here it is not useful to report the complete 
expression. However, it is interesting to report the following explicit expressions: 



W = Km£(m, V ) = SM - .log I V ~ 2{D 2 d }^ + X) ) -Irb^+f. 

(16) 



£„,(¥>) = £(!,¥>) = %>)-dlog 



n(D d ) (1 + x) 
2 2 d KipT 



2 1 + x 



log (5 



with K = —Q'{1) = 0.638..., see [20] for details. From the latter expressions one can compute the Kauzmann 
density ifK which is the solution of Y, eq (ip) = 0, and the Glass Close Packing density (fiacp which is the solution of 
£?(</?) = [20]. More generally, given a set of metastable glasses of complexity Sj, under the so-called isocomplexity 
assumption [20] , one can compute their jamming density as the solution of T,j (ip) = Sj . The results of this computation 
for binary mixtures are discussed in the main text. One should keep in mind that, as discussed above, Eqs. (16) 
incorrectly includes the mixing entropy (coming from the S(ip) term). This has to be subtracted in order to get the 
correct physical result. 



Coordination numbers 



Here we show how to compute the partial coordination numbers. The technical part of the computation follow 
closely the derivation in [20], hence the only nontrivial point is to add indices corresponding to particle types. As 
discussed in section VII. C. 3 of [20], the integral of the glass correlation function, g Miy (r), on a shell D^ v < r < 
D^u + 0(yj(fj — if) gives the number of particles of type v that are in contact with a given particle of type p for 
ip —* <Pj. A straightforward generalization of the derivation of section VII. C. 2 of [20] shows that, at the leading order 
close to contact, g^ v (r) = g^(r)[l + Q M „(r)] — g^ v [l + Q^(r)}. Then one obtains 

f D^+0( v /(A^+A„)/2) 

= QapvD^g^ / dr [1 + Q^(r)} . (17) 
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In the limit ip — > ipj, one has oc m — > and the integral of Q M1 /(r) can be easily evaluated [20]. The result is 



(18) 



where a p = lim TO ^ ^ (%Sr) • 

In the explicit case of binary mixtures, note that the equation (13) for 6 does not depend on m. Using also Eq. (14), 
we obtain 



z A a = xD^/gAAY 
zbb = D^g BB VS^ 

nd-l fl+S d 

ZAB — U AB gAB y — f 

zba = xL> AB gAB ■ 1 



(19) 



2 r 



Note that it is possible to show, using the definitions of 8 and T, that the average total coordination z = Yl^v x t* z t*u = 6, 
i.e. the packings are predicted to be isostatic irrespective of the mixture composition. 



Equation of state for liquid hard sphere mixtures 

We used a generalization of the Carnahan-Starling equation of state, which is defined by the following relation for 
the contact value of the radial distribution function [32] : 

9^) = — + (<WeM - — J {D ' d) »l , (20) 

where gpureif) is the contact value of g(r) for the pure system, given by the standard Carnahan-Starling equation [24]. 
The pressure is then given by the exact relation 

p(v) = — = 1 + j^rfH x ^ D U^ = (2i) 

Integrating this expression one obtains the entropy S(ip). The additive integration constant is fixed by the condition 
that, for f — > 0, the entropy tends to the ideal gas value S(ip) = 1 — logp — ^2 logx M . This expression includes 
the mixing entropy that must be subtracted before substituting S(<p) into Eqs. (16). 



